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The  results  of  the  efforts  for  grant  AFOSR-88-0195  are: 

1)  numerically  identified  low  dimensional  chaotic  attractors 
with  spatially  coherent  structures: 

2)  measured  the  preperties  of  the  chaos: 

3)  identified  the  sources  and  types  of  chaos; 

4)  determined  natural  coordinates  for  the  attractor  which  are 
associated  with  the  simple  spatial  patterns  in  chaos: 

5)  used  these  coordinates  to  calculate  reduced  systems  of 
equations  which  have  the  same  routes  to  chaos  and 

quali tatively-and  quatitatively-similar  strange  attractors. 
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In  the  last  seven  yezirs  there  hzis  been  a  large  group  effort  devoted  to  studying  chaos  and 
coherence  in  “near”  iategrable  pdes.  The  group  of  Alan  Bishop  of  the  Los  Alamos  National 
Laboratory,  Nick  Ercolani  of  the  University  of  Arizona,  David  McLaughlin  of  Princeton 
University,  and  Greg  Forest,  Amar  Sinha,  and  the  author  of  the  Ohio  State  University 
have  been  studying  perturbed  sine-Gordon  and  nonlinear  Schrodinger  equations.  Our 
efforts  have  been  concentrated  in  the  following  cireas: 

(1)  numerically  identifying  low  dimensional  chaotic  attractors  with  spatially  coherent 
structures; 

(2)  measuring  the  properties  of  the  chaos; 

(3)  identifying  the  sources  and  types  of  chaos; 

(4)  determining  natural  coordinates  for  the  attractor  which  zure  associated  with  the  simple 
spatial  patterns  in  chaos; 

(5)  and  using  these  coordinates  to  calculate  reduced  systems  of  equations  which  have 
the  same  routes  to  chaos  and  qualitatively  —  and  quantitatively  —  similar  strange 

■  attractors. 

This  study,  as  with  most  other  such  studies  [1,2, 3, 4, 5],  has  relied  heavily  on  numeric, 
computations  —  even  with  the  rich  structure  of  completely  integrable  systems.  However, 
because  of  this  rich  structure  we  have  been  able  to  use  numerical  experimentation  to  gain 
insight  into  the  analytical  behavior  of  the  system;  this  hais  not  generally  been  true  in  the 
pde  literature  (see,  for  example  [6, 7, 8, 9]). 

For  certain  parameter  regimes  of  the  sine-Gordon  ajid  the  nonlinear  Schrodinger  equa¬ 
tions  we  have  made  leirge  gains  in  the  first  four  areas  and  are  now  concentrating  on  the  fifth. 
From  numerical  experiments  [10]  it  is  clear  that  soliton  waveforms  are  natural  candidates, 
for  the  coordinates  of  the  attractor  and  that  only  a  very  small  number  of  these  modes  are 
sufficient.  The  soliton  perturbation  theory  of  the  1970’s  [11,12,13]  is  inadequate  to 
dress  this  calculation  of  reduced  systems  since  this  perturbation  theory  heis  great  difficu. 
in  handling  the  transitions  between  the  various  spectral  coordinates  (kink,  breather,  and 
radiation). 

Instead  we  are  exploiting  the  rich  geometric  structure  in  phase  space  of  these  integrable 
Hamiltonian  systems.  The  geometry  of  this  infinite  dimensional  phase  space  —  with  its 
separatices  and*homoclinic  orbits  —  enables  us  to  identify  precisely  the  sources  of  chans,  to 
determine  the  number  of  nonlinear  modes  required,  and  to  identify  the  transitions  between 
the  various  nonlinear  modes  [14,15].  Thus,  neeirby  integrability  provides  sufficient  control 
to  yield  very  precise  and  detailed  information  about  the  structure  of  the  chaotic  attractors 
—  much  more  detailed  information  than  can  be  anticipated  for  general  pdes. 

Of  course,  the  final  goal  of  this  lengthy  study  is  to  able  to  study  chaos  in  more  general 
systems  in  one  auid  more  space  dimensions.  For  exaimple,  our  antilytictd  cind  numerical 
results  on  homoclinic  orbits  should  be  applicable  to  a  large  cltiss  of  systems  where  this 
seems  to  be  the  oric'u  of  the  chaos.  And  our  work  on  natural  coordinates  for  the  attractor 
(i.e.,  global  and  local  bases  in  chaos)  should  be  applicable  in  systems  where  the  localized 
solutions  are  not  known  analytically.  At  present  we  are  using  what  is  known  about  chaos 
in  more  general  systems  simply  as  a  guide  to  the  areas  we  should  be  concentrating  on.  In 
the  near  future  we  expect  to  actually  apply  bur  results  to  some  of  these  systems  where  we 
will  have  to  rely  almost  totally  on  numerical  computations,  techniques,  and  results. 
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Let  me  be  specific  on  our  results  to  date  in  the  following  aresis: 

(a)  seimple  numerical  results, 

(b)  data  analysis, 

(c)  homoclinic  orbits  as  sources  of  chaos, 

(d)  geometry  of  the  integrable  sine-gordon  equation, 

(e)  reduced  system  of  odes, 

(f)  similarity  between  the  chaos  in  the  full  amd  reduced  systems,  and 

(g)  an  analytic  study  of  these  homoclinic  orbits. 

We  have  been  studying  the  driven,  damped  sine-Gordon  equation  with  periodic  bound¬ 
ary  conditions: 

Uii  —  utt  =  sinu -f- €  (aut  —  Fcoswt)  (1) 

u{x  +  L,t)  =  u(x,t),  u(x,f  =0)  =  u<„(x),  tit{x,t  =0)  =  Vin{x).  (la) 

Here  c  is  a  small  positive  parameter,  a  >  0  measures  the  strength  of  dissipation, 
the  amplitude  of  the  sinusoidal  driver  and  a>  is  its  frequency,  and  L  is  the  spatial  pe:.. 
The  initial  data  u,„(x)  is  periodic  with  period  L  and  has  one  maximum  in  this  period. 
The  initial  data  Vin(x)  is  usually  0.  (In  almost  all  the  experiments  we  have  studied  to 
date  the  initial  data  is  symmetric  and  so  the  waveform  remains  symmetric  for  all  time.) 

In  Josephson  junctions,  which  consists  of  two  superconducting  metal  layers  sopairated 
by  a  thin  insulating  oxide  layer,  which  is  small  enough  to  permit  quantum-mcchaniczd 
tunnelling  of  electrons,  different  perturbations  are  generally  used  [16,17,18]]  and  so  we  are 
now  studying 

Utr  —  utt  =  sinu  -f  6  (a«t  —  I3uxxt  —  F cosut  +  7). 

This  adds  a  dc  driver  energy  source  and  a  dissipation  which  depends  on  the  wavcnumuer. 

(a)  Sample  Numerical  Results.  We  have  been  studying  eq.  (1)  numericEilly  as  described 
in  [10,19].  By  varying  e  ,  a  ,  F  ,  u  ,  and  L  we  have  found  several  robust  attractors  with 
distinct  routes  to  chaos  (see  also  [5,20]).  The  parameter  values  we  are  concentrating  on 
are  £  =  12 ,  ea  =  0.04 ,  u  =  0.S7 ,  and  we  use  the  amplitude  of  the  sinusoidal  driver  as 
the  stress  parameter  (  0.0  <  eF  <  0.15  ).  Our  route  to  chaos  is  (see  Fig.  1  for  plots  of  these 
regions): 

cF  <  0.050  — 

spatially  fiat  -  temporally  locked  to  the  driver  (Fig.  la); 

0.050  <  eF  <  0.058  — 

spatially  one  excitation  -  temporally  locked  to  the  driver  (Figs,  lb  and  Ic); 
eF  as  0.058  — 

spatially  one  excitation  -  quasiperiodic  in  time  (Fig.  Id); 
eF  as  0.0585  — 

spatially  one  excitation  -  temporally  chaotic  (this  is  a  very  smzdl  region  which  may 
simply  be  a  result  of  our  numerical  code  —  see  Fig.  le),  the  corresponding  energy  of 
the  solution  eis  a  function  of  time  is  shown  in  Fig.  If; 

0.059  <  eF  <  0.073  — 

spatially  flat  -  temporally  locked  to  the  driver  (Fig.  Ig),  there  seem  to  be  some  very 
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small  windows  of  one  spatial  excitation  and  quasi-periodicity  in  time  which  are  not 

shown; 

0.073  <  er  — 

spatially  one  excitation  -  temporally  chaotic  (Fig.  Ih). 

Numerically  the  spatial  state  seems  to  consist  of  a  chaotic  visiting  of  three  states — a 
breather  localized  in  the  center,  a  breather  localized  in  the  wings,  2md  a  flat  state.  In  addi¬ 
tion,  the  waveform  is  low  amplitude  (~2)  and  seems  to  consist  of  two  nonlinear  modes — a 
breather  zmd  m  =  0  radiation,  and  m  =  1  and  m  =  0  radiation  (i.e.,  cos2Tmi/L  ).  We 
emphasize  that  this  chaotic  motion  occurs  for  eF  as  0.07 ,  a  value  at  which  the  system  is 
still  near  the  integrable  sine- Gordon  equation. 

This  pax2uneter  region  has  been  chosen  because  the  chaos  is  very  simple  and  the 
reduced  system  of  equations  does  a  good  job  of  explzdning  our  results.  The  details  will  be 
discussed  below  and  unless  otherwise  noted  all  discussion  and  plots  will  correspond  to  this 
pairzimeter  region. 

There  are  two  other  parameter  regions  that  aire  very  important.  The  first  is  the  sa.me 
as  the  above  except  that  L  =  24 .  The  route  to  chaos  is  similar  but  the  chaos  itself  is 
more  complicated  (10|.  That  is,  there  is  a  chaotic  bouncing  between  one  spatial  e.xcitation 
and  two  spatial  excitations  per  period.  In  addition,  if  symmetry  is  not  imposed  on  the  pde 
then  there  is  a  translational  instability  in  the  waveform. 

In  the  other  parameter  regime  there  is  a  very  unusual  period-doubling  route  to  chaos. 
This  is  where  i  =  80 ,  eoe  =  0.004 ,  and  u;  =  0.98 .  The  route  to  chaos  (as  F  increases) 
goes  from  spatially  flat  to  spatially  one  excitation  and  from  periodic  to  quasiperiodic  in 
time  [21].  However  in  this  quasiperiodicity  there  is  a  period  doubling  into  chaosl  .  That 
is,  when  we  do  a  Poincare  map  to  factor  out  the  phase  of  the  driver  then  we  see  a  perior^ 
doubling  into  chaos.  (First  there  is  just  one  point,  corresponding  to  the  periodic 
Then  there  is  a  closed  curve,  corresponding  to  the  quasi- periodic  orbit.  Then  this  do. 
curve  goes  unstable  and  becomes  a  closed  curve  which  crosses  itself,  etc.)  Then,  as  above, 
as  F  is  increased  further  only  the  flat  wave  remains  stable,  and  finally  chaos  reappears. 

It  is  our  belief  that  aJl  three  parameter  regimes  arise  from  the  same  cause  —  namely 
homoclinic  orbits.  As  will  be  described  below  we  have  verified  this  (to  our  satisfaction)  in 
the  first  region7  we  see  no  reason  to  doubt  it  in  the  second  region,  and  it  seems  reasonable 
in  the  third  region  but  much  work  remains  to  be  done. 

(b)  Data  Analysis.  In  order  to  establish  the  structure  of  these  attractors  we  have  used 
standard  diagnostics  to  analyze  our  numerical  data.  In  particulrir  the  temporal  behavior 
is  anedyzed  with  standard  diagnostics  from  dynamical  systems  theory,  such  as  time  series, 
phaise  planes,  Poincare  sections,  temporal  power  spectra,  leading  Lyapunov  exponents,  and 
correlation  dimensions.  All  temporal  diagnostics  are  consistent  with  the  above  description 
of  the  attractors. 

In  addition  we,  are  studying  the  local  behavior  of  the  orthonormal  basis  which  is  used 
to  construct  the  Lyapunov  exponents.  To  explain  this  let  us  begin  with  how  Lyapunov 
exponents  axe  calculated  both  analytically  and  numerically  [22] . 

Let  our  system,  x  =  f{x,t)  be  N  dimensional  and  consider  a  trajectory  {x(f)|t  > 
0}  where  f(0)  =  iq  .  Along  with  this  trajectory  consider  the  evolution  of  the  volume 
consisting  of  all  the  initial  conditions  that  lie  within  an  infinitesimal  ball  of  radius  € 
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centered  at  xq  .  The  evolution  of  this  ball  is  an  ellipsoid  (since  its  evolution  is  determined 
by  the  linear  part  of  the  flow)  with  principal  axes  of  length  Ci(t)  where  Ci(0)  =  «  for 
1  <  i  <  N  .  Then  the  Lyapunov  exponents,  Aj ,  are  defined  by 


Aj  =  lim  lim 

t— oo<— *0  C 


where  Ai  >  A2  >  . . . . 

Numerically,  however,  we  pick  a  time  At  eind  integrate  an  orthonormal  basis  (which 
represents  this  ball)  forward  in  time  by  At  using  the  linearized  equations.  We  then  re- 
orthonormalize  by  using  a  Gram-Schmidt  procedure  to  find  the  lengths  of  these  principcil 
axes  and  repeat  this  procedure  for  K  time  intervals.  If  we  let  the  length  of  the 
principal  axis  at  the  end  of  the  time  interval  be  6x'{k)  (where  «  =  1 )  then 


Ai 


K 


lim  -  ■—  In^xVfc) 

K-^00  KAt  ^  ^  ’ 


(3o' 


Instead  of  studying  the  Lyapunov  exponent  we  study 

ln5z‘(jt)  (3b) 

'  Jbsico 

directly  as  a  fxmction  of  {k  —  kQ)At .  (We  call  these  the  local  Lyapunov  exponents  —  ILe’s  '' 
The  reaison  for  our  choice  can  be  seen  in  Figs.  2.  In  Fig.  2a  we  show  the  convergen- 
(In5i^(fc))/At  as  a  function  of  kAt .  There  is  one  positive  Lyapunov  exponent,  two  zero 
Lyapxmov  exponents  (if  symmetry  is  used  in  the  code  then  there  is  only  one),  and  the 
other  five  Lyapunov  exponents  (that  were  calculated)  have  the  seime  negative  value.  In 
Fig.  2b  we  plot  ^^ln5i^(fc)  as  a  function  of  {k  —  ^o)At  where  the  queisiperiodic-like 
character  of  the  chaos  can  be  eeisily  seen.  Note  that  this  repeated  signature  —  of  length 
~  50  corresponds  to  the  random  fluctuations  of  the  breather-like  waveform  (as  in  Fig. 
Ih)  as  the  waveform  dies  out  aind  then  recovers  either  as  before  or  shifted  by  L/2  .  Thus 
the  average  length  of  these  signatures  gives  the  time  for  this  recovery  —  although  it  does 
not  indicate  whether  the  new  breather-like  waveform  is  shifted  with  respect  to  the  old  or 
not.  In  addition,  the  occasional  huge  jumps  in  this  local  Lyapunov  exponent  are  caused 
by  the  flat  waveform  persisting  for  a  much  longer  period  of  time.  In  the  ode  systems  this 
corresponds  to  the  trajectory  circling  around  an  unstable  fixed  point  for  a  few  rotations 
before  zooming  off  again.  In  Fig.  2c  we  again  plot  a  function  of  (fc  —  kQ)At 

but  for  the  6  dimensional  sine- Gordon  ode.  At  least  to  the  naked  eye  these  curves  are  very 
similar  —  even  to  having  the  same  huge  jumps. 

There  eire  many  advantages  to  studying  these  local  Lyapunov  exponents: 

1)  These  ILe’s  are  functions  of  time  and  so  cm  exhibit  much  interesting  behavior  over 
the  whole  trajectory  Jis  opposed  to  the  Lyapunov  exponents  which  are  a  time  average 
over  the  whole  trajectory. 
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2)  These  ILe’s  axe  ordered  so  that  the  leading  one  or  ones  show  the  exponential  rate  of 

separation  which  is  the  characteristic  of  cheu)s.  In  fact  for  the  sine-Gordon  pde  and  for 
our  particular  parzoneter  regime  Ai  >  0 ,  =  A3  =  0  ,  and  A4  =  A5  =  . . .  =  Ag  <  0 

(which  is  all  the  Lyapunov  exponents  we  have  calculated  since  each  one  takes  as  long 
to  solve  as  the  original  pde).  The  first  ILe  (shown  in  Fig.  2b)  and  the  second  and 
third  ones  have  the  same  characteristic  signature  amd  look  quite  similar  (although 
their  slopes  are  quite  different).  The  remaining  ILe’s  all  look  similar  to  the  ILe’s  for 
non-chaiotic  waveforms. 

3)  These  ILe’s  give  am  overall  picture  of  the  trajectory  as  opposed,  for  example,  to  plots 
of  the  j***  component  of  the  trajectory. 

4)  We  can  study  statisticad  properties  of  these  ILe’s  amd  so  compare  chamtic  trajecto¬ 
ries  between  the  sine-Gordon  pde,  the  nonlineau:  Schrodinger  pde,  the  sine-Gordon 
odes,  amd  the  nonlinear  Schrodinger  odes.  This  includes  studying  correlations,  power 
spectra,  etc.  in  order  to  determine  quamtitative  measures  of  the  differences  between 
different  trajectories. 

The  spatial  behavior  is  anailyzed  with  the  inverse  scattering  transform  which  measures 
quantitatively  both  the  soliton  and  radiation  nonlinear  modes  that  comprise  the  attractor 
at  a  given  time  [19] .  With  L  =  12  we  find  that  the  number  of  appreciably  excited  nonlin¬ 
ear  modes  is  2  at  the  onset  of  chaos  amd  increases  slowly  with  increasing  cF  .  The  chamtic 
attractor  itself  can  be  described  by  these  modes  imdergoing  collision,  amnihilation,  nucle- 
ation,  amd  energy  tramsfer  between  coherent  (locadized)  states  and  radiation  (extended) 
states.  These  transitions  involve  what  we  call  the  “crossing  of  homocUnic  states”  and  are 
shown  in  Figs.  3.  Notice  how  the  spectrum  (the  wiggly  line)  oscillates  between  a  cross  (for 
exaimple  at  f  =  25,036 )  amd  a  gap  (for  example  at  f  =  25,048 ).  The  homoclinic  orbit 
corresponds  to  the  spectral  configuration  between  these  two  times  where  the  gap  closes 
and  the  two  pieces  of  spectrum  just  touch.  As  time  evolves  these  two  pieces  of  spectr 
merge  into  one  and  push  the  spectrum  out  (approximately)  perpendicular  to  the  ci. 
indicated  by  the  dots  (this  is  discussed  in  detail  in  [27]).  The  gap  configuration  always 
corresponds  to  higher  energy  and  to  a  more  focused  waveform  (i.e.,  more  localized). 

(c)  Homoclinic  Orbits  as  Sources  of  Chaos.  Often  the  origin  of  chaos,  amd  of  its 
associated  sensitivity,  cam  be  traced  to  homoclinic  orbits  of  a  nearby  integrable  system. 
This  source  of  chaos  is  well  understood  in  odes.  For  example,  in  the  damped,  driven 
pendulum  the  unperturbed,  infinite  period  separatrix  is  a  homoclinic  orbit.  The  motion 
of  the  pendulum  neau'  its  inverted  position  is  extremely  sensitive  to  the  sinusoidal  driver. 
Thus,  smadl  perturbations  cam  cause  the  pendulum  to  either  return  to  its  original  potential 
well  or  to  be  kicked  into  the  next  potential  well.  Ramdom  and  chaotic  behavior  results 
as  a  consequence  of  this  sensitivity  [23|.  In  the  pde  caise  detailed  geometric  information 
about  infinite  dimensional  phaise  space  is  generadly  not  available.  However,  for  the  special 
integrable  pdes  of  soliton  theory  such  geometric  information  is  becoming  available.  For 
the  Korteweg-de  Vries  equation  under  periodic  boundary  conditions  a  rather  complete 
geometric  theory  exists  [24].  It  is  interesting  to  note  that  this  periodic  KdV  system  hats  no 
homoclinic  orbits,  no  instabilities,  and  no  separatices — amd  requires  very  strong  stresses  to 
become  chaotic.  On  the  other  hand  the  periodic  sine-Gordon  equation  possesses  a  much 
richer  geometric  structure  [25]  than  KdV.  The  sine-Gordon  equation  has  many  homoclinic 
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states  throughout  its  phase  space  whichprovide  a  source  of  chaos  when  the  system  is  driven; 
and,  indeed,  the  chaotic  structure  is  very  rich. 

Homoclinic  orbits  can  occur  whenever  the  spectrum  has  a  critical  point  [15].  That  is, 
the  spectrum  is  all  the  values  of  A  where  an  amalytical  function  called  the  discriminant, 
A(A),  is  real  and  — 2  <  A  <  +2 .  A  critical  point  Ac  is  wherever  A(Ac)  =  ±2  and 
A'(Ac)  =  0 .  If  a  critical  point  satisfies  3?(Ac)  >  0  and  3f(Ac)  >  0  then  it  can  have 
exponential  growth  [15].  We  have  explicitly  calculated  the  homoclinic  orbit  which  seems 
to  be  the  source  of  the  chaos.  It  is  shown  in  Fig.  4a  and  a  representative  plot  of  the 
chaotic  trajectory  from  Fig.  Ih  is  shown  in  Fig.  4b.  As  a  more  accurate  comparison  we 
plot  u(x  =  0,  t)  for  both  the  homoclinic  orbit  and  the  chaotic  trajectory  in  Fig.  4c. 

Thus,  the  nonlinear  trzuisform  has  established  that  imperturbed  homoclinic  states, 
which  2U'e  potential  sources  of  sensitivity,  are  indeed  present.  We  emphasize  that  this 
explicit  and  precise  identification  of  the  homoclinic  states  and  their  crossings  is  only  pos¬ 
sible  because  we  are  in  a  near  integrable  situation  and  have  a  precise  tool  (the  non’ 
transform)  for  their  detection. 

(d)  Geometry  of  the  Integrable  Sine-Gordon  Equation  This  theoreticaJ  work  [14,15] 
may  be  summarized  geometrically  as  an  extension  of  the  2-dimensional  phase  portrait  for 
the  pendulum  into  the  infinite  dimensional  phase  space  of  the  integrable  periodic  sine- 
Gordon  equation,  and  physically  as  an  extension  of  the  classical  description  of  a  long 
wavelength,  linear  instability  of  trivial  solutions  to  both  a  linear  and  a  global  description  of 
similar  instabilities  for  arbitrary  N-phase  solutions.  It  is  this  basic  theory  that  has  enabled 
us  to  detect — with  the  nonlinear  spectral  measuring  device — the  presence  of  homoclir’ 
crossings  in  the  numerical  data  for  the  chaotic  attractor. 

The  theoretical  understanding  we  have  obtained  about  the  Benjamin-Feir  unstable 
modes  has  given  an  analytical  underpinning  to  the  numericail  calculations  of  the  presence 
of  homoclinic  orbits  in  the  chaotic  waveforms.  Additionzdly  we  have  been  able  to  explicitly 
write  down  the  homoclinic  orbit  which  seems  to  be  the  cause  of  the  chaos  in  the  parameter 
region  with  L  =  12  as  described  above.  (This  solution  is  no  easy  feat  since  it  involves  the 
Jacobi  eta,  theta,  and  zeta  functions!) 

(e)  Reduced  System  of  Odes.  Since  the  geometry  of  the  sine-Gordon  phase  space  is 
quite  complicated  [14,15],  we  have  constructed  reduced  systems  of  odes  to  approximate 
the  infinite  dimensioned  pde  using  both  Fourier  modes  emd  Lumley  modes  [26].  (Lumley 
modes  au-e  orthogonal,  coherent  modes  in  chaos  which  axe  “best”  in  some  least-squeires 
sense).  Preliminary  work  has  begun  by  Forest,  McLaughlin,  and  Sinha  on  e.xpanding  the 
waveform  in  a  nonlinear  basis  (using  the  geometry  of  the  sine-Gordon  phase  space).  Early 
results  axe  promising.  However  our  work  will  continue  to  focus  on  linear  bases. 

Three  reduced  systems  have  been  studied  in  great  detail.  The  first  is  a  reduced  system 
of  ordinary  differential  equations  directly  taken  from  the  sine-Gordon  pde.  That  is,  we 
assume 

OO 

u(i,f)  =  ^njk(0efc(x)  .  (4) 

k=Q 
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(5) 


and  truncate  this  linear  combination  to  solve  the  2iV***  order  sine-Gordon  ode  system 


d*  =  ak{e'l,tk)  —  ysin(^  ak6k),^kj  —  €aak  +  €rsina;i  (l,efe) 

for  =  0, 1, 2, . . . ,  iV  —  1 .  There  are  two  possibilities  here:  1)  leave  the  problem  as  it  is  and 
calculate  the  inner  product  containing  the  sine  term  at  each  time  step;  or  2)  expand  the 
sine  term  into  the  first  two  or  three  terms  of  its  Taylor  series  so  that  this  inner  product  can 

be  calcTilated  once  and  for  all.  The  second  approach  has  the  virtue  of  saving  a  significant 

amount  of  computer  time.  The  first  approach  has  the  advantage  of  being  more  accurate 
and  allowing  any  number  of  basis  functions  to  see  if  some  sort  of  convergence  is  obtained 
as  JV  '■+  oo .  We  are  using  both  approaches.  We  are  also  using  two  different  choices  for 
the  et’s ,  namely,  cosine  modes  and  Lumley  modes. 

.The  second  is  to  derive  a  perturbed  nonlinear  Schrodinger  envelope  equation  from  e 
(1).  Letting 

w  =  1  —  ew,  ea  =  2u>a,  eF  =  8V^  u>f ,  L  =  v^2«I>X, 

X  =  n/2^x,  T  =  eut,  u(x,  t)  =  2v^  [B{X,  T)  exp“^‘  +  c.c.]  +  0(e) 

and  keeping  the  first  two  terms  in  the  Taylor  series  expansion  for  sin  u  ,  then  the  slowly 
varying  envelope  B(X,  T)  satisfies 

—iBr  +  Bxx  +  (1-Sl^  —  \)B  —  ictB  —  f.  (7^ 

We  have  achieved  two  things  by  this  reduction.  First,  we  have  preserved  the  pertur^ 
integrable  structure  since  the  unperturbed  pde  of  eq.  (7)  (i.e.,  d  =  f  =  0  )  is  the  integrable 
nonlinear  Schrodinger  equation.  Second,  we  have  factored  out  the  frequency  of  the  dhiver, 
u; .  Thus,  periodic  solutions  of  the  sine-Gordon  equation  correspond  to  fixed  points  of 
the  nonlinear  Schrodinger  equation,  quasiperiodic  solutions  of  the  sine-Gordon  equation 
correspond  to  periodic  solutions,  which  are  incommensurate  with  u  ,  of  the  nonlirear 
Schrodinger  equation,  and  chaotic  solutions  of  the  sine-Gordon  equation  correspond  to 
chaotic  solutions  of  the  nonlinear  Schrodinger  equation. 

We  now  maike  a  further  approximation  eind  truncation  based  on  the  numerically  ob¬ 
served  low  number  of  nonlinear  modes,  namely 

N-\ 

B{XyT)  =  ^kit)  cos  kX  (8) 


(where  the  i4’3  are  complex).  We  then  obteun  an  explicit  system  of  ordinary  differential 
equations  in  Ak  using  Mathematica.  At  present  we  have  used  N  <  5  . 

In  particular,  tor  N  =  2  we  obtain  the  4^**— order  system  of  odes 

— zcr  +  ((c(^  + -|6(^  —  l)c -f- -(c6 -f- c6)6  =  idc  —  f, 

3  “  -  -  . 

—  ibx  -|-  (|c)^  -I-  —  (1  -f  -f-  (cb  -I-  cb)c  =  iab 
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which  is,  in  fact,  an  integrable  Hamiltonian  system  [27] .  Its  constamts  of  the  motion  are 

13  1  1  (^0) 

H  =  i|cr  +  l^ncl^  +  ~\b\<  -  ~(1  +  -  |cp  +  +  b^c^). 

The  fixed  points  of  the  ode  correspond  qualitatively  to  the  periodic  solutions  of  the  pde 
before  the  onset  of  chaos;  and,  in  chaos,  the  three  unstable  fixed  points  of  the  ode  which 
seem  to  be  the  basis  for  the  attractor  correspond  in  the  pde  to  a  breather  localized  in  the 
center,  a  breather  localized  in  the  wings,  and  a  flat  state  (see  Fig.  Ih).  In  addition,  the 
homoclinic  crossings  of  the  pde  correspond  quite  accurately  to  the  zeroes  of 

h  =  (11) 


Although  this  ode  tnmcation  cannot  be  expected  to  yield  quantitative  agreement 
with  the  pde  it  does  model  most  of  the  apparent  features  of  the  route  into  clxaos  and 
the  behavior  in  chaos  of  the  pde.  This  is  somewhat  surprising  since  the  approximations 
necessary  to  obtain  eq.  (7)  require  replacing  sinu  by  u  —  u^/6  (even  though  max{u}  >2) 
and  neglecting  the  next  term  in  the  expansion  in  e  ,  namely  0{y/e) . 

(f)  Similarity  Between  the  Chaos  in  the  Full  and  Reduced  Systems.  To  verify  that 
the  chaos  in  the  full  sine-Gordon  pde  and  in  various  lower  order  truncations  axe  similar  we 
need  to  show  that  the  strange  attractor  is  similat  in  some  sense.  As  discussed  abc 
method  to  determine  the  similarity  between  two  chaotic  trajectories  is  to  compme 
local  Lyapunov  exponents.  This  comparison  is  shown  in  Figs.  2b  and  2c.  To  make  i.... 
comparison  more  qualitative  we  have  taken  the  power  spectra  of  these  two  figures  and  plot 
them  in  Figs.  5.  We  are  presently  working  with  Mark  Berliner  of  the  Statistics  Department 
in  trying  to  determine  how  precisely  we  can  say  that  two  chaotic  trajectories  are  “similar’’. 
This  includes  studying  correlations,  power  spectra,  etc.  in  order  to  determine  quantitative 
measures  of  the  differences  between  different  trajectors.  Our  mmn  emphasis  is  using  the 
ILe’s  since  these  give  a  basis  free  view  of  the  strange  attractor.  At  present  we  can  say 
that  certain  trajectories  are  “closer”  tham  others  but  we  still  need  to  develop  a  technique 
for  calculating  error  bounds  so  that  we  can  show  that  two  trajectories,  corresponding  to 
the  same  penameters  but  different  initial  conditions,  are  “very  similar”.  In  this  way  we 
could  check  whether,  for  example,  chaotic  trajectories  in  the  sine-Gordon  odes  converge  to 
chaotic  trajectories  in  the  sine-Gordon  pde  as  we  increase  their  dimension.  {They  certainly 
seem  to  be  “close”.)  In  this  way  we  could  show  (at  least  numerically)  that  the  geometry 
underlying  the  sine-Gordon  odes  is  the  cause  of  the  chaos  in  the  sine-Gordon  pde. 

As  a  second  treasure  of  the  similarities  between  the  various  systems  we  are  using 
the  bifurcation  sequences  of  the  various  ode  truncations  which  we  are  calculating  using 
a  public  domain  bifurcation  pacliage  AUTO  [28]  as  well  as  a  bifurcation  code  written  by 
Xiong  to  provide  backup  support  for  AUTO  for  verification,  to  use  in  place  of  AUTO  when 
AUTO  seems  not  to  be  working  correctly,  and  to  study  in  more  detail  the  structure  of  this 
bifurcation  sequence.  In  Figs.  6a  and  6b  we  compare  the  bifurcation  sequence  for  the 
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6‘**— order  nonlinezir  Scorodinger  ode  and  the  6‘**— order  sine-Gordon  ode  respectively. 
We  cannot  follow  the  Hopf  bifurcation  in  the  sine-Gordon  ode  since  it  corresponds  to  a 
quasi- periodic  orbit,  but  otherwise  the  two  sequences  are  very  similar.  (The  norms  (the 
vertical  axes)  are  not  the  same  in  both  plots  and  so  the  curves  look  different  —  but  the 
criticad  points  all  occur  at  approximately  the  same  values  of  eF ).  In  particular,  as  the 
order  of  the  sine-Gordon  ode  increases  the  critical  points  seem  to  approach  the  boundaries 
of  the  sine-Gordon  pde.  Even  for  the  6***— order  ode  the  accuracy  is  to  within  ~  10%  for 
all  the  critical  points  shown.  In  pzirticular,  whereas  the  4***— order  ode  has  a  subcritical 
Hopf  bifurcation  (and  so  does  the  4‘^— order  nls  ode)  cdl  higher  order  systems  have  a 
supercritical  Hopf  bifurcation  which  corresponds  to  the  quasi-periodicity  in  the  pde  at 
cF  ~  0.058 .  Also,  the  extent  (in  F  )  of  these  quasi-periodic  solution  agrees  very  well  with 
the  range  found  numerically  in  the  pde.  Thus  we  axe  fully  confident  that  the  geometry 
imderlying  the  ode  systems  can  explain  the  transition  to  chaos  in  the  pde. 

(g)  An  analytic  study  of  these  homoclinic  orbits.  Dave  McLaughlin  and  Steve 
Wiggins  are  working  with  Chuyu  Xiong  and  me  in  attempting  to  use  singular  perturbation 
theory  to  understand  the  nature  of  the  homoclinic  orbits  which  are  the  sources  of  the  chaos. 
This  work  is  still  preliminary  but  Xiong  and  I  are  developing  better  numerical  techniques  to 
study  periodic  orbits  which  axe  very  close  to  the  homoclinic  orbit.  (We  have  not  been  able 
to  develop  any  techniques  to  solve  for  the  homoclinic  orbit  directly  since  the  fixed  point 
corresponding  to  this  homoclinic  orbit  is  2-dimensionally  unstable.)  With  our  numerical 
results  they  are  able  to  determine  how  to  do  their  perturbation  expansions. 

Publications 

Two  papers  have  misen  so  fax  from  this  study.  In  the  first  [29]  we  did  an  analytical 
and  numericaJ  perturbation  analysis  on  the  nonlineax  Schrodinger  equation  for  ea  =  0.04  , 
L  =  12  cind  L  =  24 ,  and  for  various  uj  and  F  .  This  entailed  calculating  the  waveform  in  a 
perturbation  expansion  and  using  this  in  stability  calculations.  This  allowed  us  to  validate 
the  solutions  we  were  seeing  numerically  in  the  pde  and  to  obtain  our  first  bifurcation 
sequence  (similar  to  Figs.  6)  for  comparison.  It  also-allowed  us  to  understand  how  the 
bifurcation  sequence  changes  as  a  and  u  vary.  However  this  paper  was  concerned  strictly 
with  the  physical  plane  and  never  made  use  of  the  spectrail  transformation  underlying  the 
sine-Gordon  pde  (i.e.,  the  underlying  geometry  of  the  pde). 

In  the  second  paper  [27]  the  first  attempt  'vaa  made  to  study  the  geometry  of  the 
pde  explicitly.  In  this  study  we  attempted  to  understand  the  transition  to  chaos  for 
the  parameter  regime  L  =  12  .  ea  =  0.04 ,  eind  u)  =  0.87  (although  the  analysis  is 
certainly  applicable  for  other  regimes).  This  was  done  by  studying  the  4^**— order  nonlinem 
Schrodinger  equation,  eq.  (9)  ode  both  numerically  and  analytically.  The  unperturbed 
system  is  completely  intv^grable  and  so  has  a  well-understood  geometric  structure.  We 
were  able  to  show  that  the  stable  fixed  points  of  the  perturbed  system  corresponded  to 
the  flat  and  periodic  solutions  of  the  sine-Gordon  pde.  We  were  also  able  to  show  that 
the  chaos  in  the  ode  arose  from  the  fixed  points  corresponding  to  these  periodic  solutions 
becoming  unstable.  Finally  we  could  identify  a  homoclinic  orbit  in  the  sine-Gordon  pde 
with  a  function  of  the  two  constzuits  of  motion  in  the  ode  system  and  thus  construct  the 
homoclinic  orbit  we  believed  to  be  the  source  of  the  chaos  explicitly  for  the  ode.  We  could 
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then  show  that  this  ode  homoclinic  orbit  was  crossed  repeatedly  in  chaos.  Of  course,  this 
study  was  carried  out  in  a  very  low  dimensional  ode  and  so  it  was  certainly  not  obvious 
that  the  chaos  in  the  pde  could  be  understood  by  this  simplified  model  (and  we  made  that 
point  in  the  paper). 

At  present  we  are  writing  up  our  results  in  comparing  the  sine-Gordon  pde  with  all 
of  the  various  truncations  discussed  above.  This  will  leave  no  doubt  that,  in  fact,  the 
geometry  of  the  low  dimensional  ode  truncations  (as  long  as  they  are  at  least  6^**— order  ) 
is  sufficient  to  explain  the  chaos  in  the  pde  by  showing  the  convergence  to  the  limit  of 
the  sine-Gordon  odes  to  the  sine-Gordon  pde  [30].  Additionally  it  will  demonstrate  new 
methods  for  showing  the  “closeness”  between  various  strange  attractors  which  will  allow 
us  to  quickly  determine  the  similarity  between  the  chciotic  responses  of  various  systems. 

Finally  we  are  writing  up  our  work  on  the  explicit  calculation  of  the  homoclinic  orbit 
which  seems  to  be  the  cause  of  the  chaos  in  the  pde  [31|.  This  will  include  a  discussion  of  the 
various  solutions  that  arise  from  one  and  two  Backlund  transformations  of  the  McCumber 
solution  (i.e.,  spatially  flat)  and  their  stability.  Both  breather-like  (i.e.,  symmetric)  and 
kink-like  (i.e.,  asymmetric)  waveforms  will  be  studied. 

Invited  Talks 

Invited  talks  on  these  subjects  were  given  at: 

•  Conference  on  Nonlinear  Transport  Properties  in  Condensed  Matter  Physics,  July, 
1989,  at  the  Technical  University  of  Denmark. 

•  IMACS  International  Conference  on  Computational  Physics,  June  1990,  at  the  Uni¬ 
versity  of  Colorado. 
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Figure  Captions 


1.  Numerically  computed  solutions  of  the  perturbed  sine-Gordon  equation  for  L  =  12 , 
CQ  =  0.04 ,  (jj  =  0.87  and  for  various  values  of  cF .  (a)  eF  =  0.48  ,  (b)  eF  =  0.50  ,  (c) 
eF  =  0.57 ,  (d)  eF  =  0.58  ,  (e)  eF  =  0.586 ,  (f)  the  energy  as  a  f\inction  of  time  for  the 
waveform  shown  in  (e)  to  make  the  evolution  clearer  (also  shown  axe  the  averages 
in  I  of  u  and  uj  as  a  function  of  time),  (g)  eF  =  0.72  ,  and  (h)  eF  =  0.75  . 

2.  Plots  of  the  leading  Lyapunov  exponents  for  the  run  with  parameters  L  =  12  ,  ea  = 
0.04 ,  eF  =  0.075 ,  and  u  =  0.87 .  (a)  The  Lyapunov  exponents  are  the  limiting 
values  of  these  curves  (i.e.,  as  <  — ♦  oo  ).  (b)  The  local  Lyapunov  exponent  for  the  the 
sine-Gordon  pde.  (c)  The  local  Lyapunov  exponent  for  the  6  dimensional  sine-Gordon 
ode. 

3.  Spectra  for  the  run  with  parameters  L  =  12,  cot  =  0.04 ,  eF  =  0.075  ,  emd  oj  =  0.87  . 

■  Some  crosses  and  gapes  are  noted  in  the  plots,  (a)  The  sine-Gordon  spectra  at  various 

times,  (b)  The  plot  of  u{x,t)  as  a  function  of  t  corresponding  to  the  times  of  the 
spectra. 

4.  (a)  The  homoclinic  orbit  which  appears  to  be  the  source  of  the  chaos  for  the  run 
shown  in  Fig.  Ih.  (b)  An  enlarged  view  of  the  chantic  trajectory  to  the  same  scale  as 
(a),  (c)  A  comparison  of  the  homoclinic  orbit  with  the  chaotic  trajectory  by  plotting 
u(x  =  0,  t)  for  each  (the  solid  line  is  the  homoclinic  orbit  and  the  dashed  line  is  the 
chaotic  trajectory). 

5.  Power  spectra  of  the  leading  local  Lyapunov  exponent  corresponding  to  the  same  runs 
as  in  Fig.  2.  (a)  The  pde;  (b)  the  ode. 

6.  Bifurcation  sequence  for  the  order  nonlinear  Schrodinger  ode  (a)  and  the 
6'**— order  sine-Gordon  ode  (b).  (The  norms,  i.e.,  the  vertical  axes,  axe  not  the  sar¬ 
in  both  plots  and  so  the  curves  look  different  —  but  the  critical  points  all  occur  a 
approximately  the  same  values  of  eF .  Also  the  horizontal  axes  axe  not  to  the  same 
scale) 

7.  The  time  evolution  (a)  and  the  corresponding  spectrum  (b)  for  £  =  80  ,  ea  =  0.004  , 
eF  =  0.0043,  and  u;  =  0.87 . 

8.  (a)  The  time  evolution  of  a  solution  to  the  exact  nonlinear  Schrodinger  pde  using  24 
points,  (b)  The  time  evolution  of  the  spectrum  during  the  time  that  the  solution  begins 
to  move.  At  each  time  the  top  two  plots  show  the  waveform  as  its  magnitude,  phase, 
real  part,  and  imaginary  part.  The  main  plot  shows  the  spectrum.  The  spectrum 
is  the  solid  lines  along  with  the  horizontal  axis  (’G’  means  the  spectrum  is  in  a  gap 
configuration,  ’C’  means  it  is  in  a  cross  configuration,  and  ’S’  rneems  the  spectrum  has 
split  aisymmetrically).  The  small  crosses  axe  where  the  discriminant  is  real  but  less 
than  —2  and  the  small  pluses  are  where  the  discriminant  is  real  but  greather  than 
-f-2 .  The  bottom  plot  shows  the  values  of  the  discriminant  on  the  real  axis. 
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